Phase space reduction of the one-dimensional Fokker-Planck (Kramers) equation 
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A pointlike particle of finite mass m, moving in a one-dimensional viscous environment and 
biased by a spatially dependent force, is considered. We present a rigorous mapping of the Fokker- 
Planck equation, which determines evolution of the particle density in phase space, onto the spatial 
coordinate x. The result is the Smoluchowski equation, valid in the overdamped limit, m — ► 0, with 
a series of corrections expanded in powers of m. They are determined unambiguously within the 
recurrence mapping procedure. The method and the results are interpreted on the simplest model 
with no field and on the damped harmonic oscillator. 

PACS numbers: 05.40.Jc, 87.10.Ed 



I. INTRODUCTION 

The Brownian motion of a particle in a confined sys- 
tem represents an essential model used in description of 
stochastic transport through quasi one-dimensional (ID) 
systems, e.g. channels in nanomaterials, pores or fibers 
in biological structures. In a ID system, the trajectory 
x(t) of a Brownian particle is described by the Langevin 
equation 



rax + jx 4- d x U (x) = f(t) 



(1.1) 



Here, m denotes mass of the particle, U(x) represents 
the driving potential, 7 is an effective friction coefficient 
and f(t) is the stochastic force, satisfying the usual con- 
ditions on averaged values, (/(<)) = and (f(t)f(t')) = 
2jkBTS(t — t')\ T is the temperature and ks the Boltz- 
mann constant. The corresponding phase space density 
p(x, v, t) of the particle satisfies the Fokker-Planck (FP), 
or Kramers (kinetic) equation, 

dt + vd x a v p(x, v, t) 

m J 

= -^d v e-P mv2 / 2 d v e? mv2 / 2 p(x,v,t), (1.2) 

where x is the spatial coordinate, v denotes its velocity 
and p = 1/ksT is the inverse temperature. 

Solutions of the Langevin equation, as well as the cor- 
responding kinetic equation, are studied over almost a 
century [l|-|3|. Still, motion of a particle in confined ge- 
ometries represents usually a complicated problem, re- 
quiring next reductions of the used description. Due 
to simplicity, and often for relevance, mainly the over- 
damped limit is studied. The mass dependent term in Eq. 
([l.ip . rax, is then considered negligible and the particle's 
spatial density p(x, t) is governed by the Smoluchowski 
equation, 

dtp(x, t) = D d x e-^d x e^p(x, t) , (1.3) 



containing no information about the mass of the parti- 
cle; Do denotes the diffusion constant. Instead of the full 
phase space, one works with only the spatial coordinate 
x. Of course, the question is the price for this simpli- 
fication, as well as possibility of also including properly 
inertia of the massive particles in the reduced (real space 
only) description of the Brownian dynamics. Recent 
studies showed its importance in understanding rectifi- 
cation of the transport in ratchets H, or its influence on 
the effective diffusion coefficient [5j-|7| in a narrow chan- 
nel. For our purpose to demonstrate the phase space re- 
duction, we will deal here with only the ID FP (Kramers) 
equation (|1.2I) . 

The eq. (I1.2j) is reducible to the Smoluchowski eq. 
(II. 3|) in the limit m — > 0. The reduction procedure [8|- 
|l0| is based on an instant thermalization of the particle's 
velocity after any move in the x direction in the case of an 
infinitely small mass. The situation resembles derivation 
of the Fick- Jacobs (FJ) equation [ill [l^. 



d t p(x,t) = d x A(x)d, 



' A{x) • 



(1.4) 



appearing as the result of the dimensional reduction of 
the diffusion equation in a 2D channel with varying cross 
section A(x), onto the longitudinal coordinate x, if equi- 
libration in the transverse (y) direction is instant. The 
function p(x,t) denotes the linear (ID) density of the 
particle. Of course, as in reduction of diffusion to the FJ 
equation, that of (|1.2[) to (|1.3|) as m is a singular 
limit, and must be handled with care, but from the view- 
point of the reduction of flTTTJ) , no such caveat is needed. 

Recently, an exact mapping procedure was proposed 
(l3l [Til , enabling us to also derive the corrections to the 
FJ equation (jl.4p . which are necessary, if the transverse 
equilibration is not instant. The key was to introduce 
anisotropy of the diffusion constant in the diffusion equa- 
tion, 

d t p(x,y,t) = \ dl + -dy) p(x,y,t), (1.5) 
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governing the 2D (spatial) density p(x,y,t). For e — » 0, 
the infinitely fast transverse diffusion immediately flat- 
tens the y profile of p(x,y,t) — > p(x,t) — p(x,t)/A(x). 
Then integration of Eq. (|1.5p over the cross section, to- 
gether with the reflecting boundary conditions satisfied 
at the hard walls, results in the FJ eq. (|1.4[) . In the 
case of a slower transverse diffusion, the mapping proce- 
dure generates a series of corrections to the FJ equation 
controlled by e > 0, developing the y profile of p(x, y, t), 
which is already curved depending on the flux and geom- 
etry of the channel. 

The procedure was extended to mapping of diffusion 
in an external field U(x, y), e.g. for diffusion in a channel 
with soft walls 15|, where the particle is kept near the 
x axis by the parabolic potential, (3U{x,y) = a(x)y 2 ; 
a(x) represents the varying stiffness of the walls. The 
equation to be mapped onto the x coordinate is the 2D 
Smoluchowski equation, 

(d t - d x e~ a ^y 2 d x e a ^y 2 ) P (x, y> t) = 

= -dye-^^ dye 01 ^ p(x, y t t). (1.6) 

Integrating over y and applying the mapping scheme 
gives the mapped ID equation in an extended Smolu- 
chowski form, 

d tP (x,t) = d x e~ v W \l + e Z{x,d x )\ d x e v ^p(x,t), 

(1.7) 

governing the ID density p(x,t), where V{x) stands for 
an effective potential and Z is the correction operator 
sought as an expansion in the small parameter e. The 
potential V(x) and the coefficients of Z depend on a{x) 
and both are fixed unambiguously within the recurrence 
mapping procedure. 

The central idea of this paper is a conjecture that 
the FP eq. (| 1 . 2|) can be reduced to an extended 
Smoluchowski-like ID form, governing the spatial den- 
sity p(x, t), after integration over v and applying a similar 
mapping scheme. The velocity v thus represents a " trans- 
verse" coordinate instead of y, with the mass to playing 
the role of the small parameter e. Then the Smoluchowski 
equation (11.31) should be obtained in the limit to — >• 
[H, E|. Performing the recurrence procedure, a series of 
corrections to this equation in powers of to would be de- 
rived. Then the final mapped equation would also respect 
inertia of the Brownian particles, although working only 
in real space. 

Use of the mapping procedure developed for diffusion 
is not straightforward; the left-hand side oper- 
ator of the FP equation (|1.2p has a different structure 
than that of the diffusion (ll.5[) or Smoluchowski equation 
(11.6[) . Still, there is a way to apply the general scheme 
of the mapping in this case and to perform reduction of 
the phase space onto the real space in the way described 
above. Presentation of this algorithm is the primary aim 
of this study. 



The result of the mapping of the FP equation (|1.2I) is 
again an equation of the form (|1.7[) with V(x) = (3U{x) 
and e replaced by to. In the limit of stationary flow, 
this equation can be simplified by subsequent reduction 
of the operator 1 — mZ to a function D(x), a spatially 
dependent effective diffusion coefficient [12| ll7| , 

d t p(x, t) = d x e-^D(x)d x e^p(x, t) . (1.8) 

The leading term of D(x) is proportional to mU"{x) and 
the whole series of corrections to the Smoluchowski equa- 
tion can be summed up, giving 



D(x) = D Q 



1 - y/l -4mU"(x)/-f 2 
2mU"(x)/-/ 2 



(1.9) 



with Dq = 1/7/3, if the higher derivatives of U(x) are 
neglected. 

In the following section, we analyze how to apply the 
mapping scheme for reduction of the phase space in the 
FP equation (|1.2j) . In Sect. Ill, our considerations are 
verified on an exactly solvable model, the FP equation 
with no field, U(x) — 0. Analysis of this example helps 
us to construct the recurrence scheme for calculation of 
the series of corrections to the zeroth order approxima- 
tion, the Smoluchowski eq. (|1.3I) . in the small parame- 
ter to. Finally, the complete mapping procedure for an 
arbitrary potential U(x) is presented in Sect. IV. The 
mapped equation of type (11.7[) . as well as the formula 
(11.9[) is derived, and checked on the damped linear har- 
monic oscillator. 



II. PRELIMINARY CONSIDERATIONS 

The key points of the mapping procedure, as formu- 
lated for diffusion 13, 14|, are recalled in this Section. 
Based on physical considerations, we adjust the general 
scheme of the mapping for dimensional reduction of the 
FP equation (|1.2j) to this situation. 

The mapping procedure represents a consistent tran- 
sition from a fine grain to a coarse grain description of 
some evolution process. The process is described in de- 
tails by some partial differential equation (PDE), gov- 
erning the density of particles p(r,t) in the full space, 
defined by the coordinates r. The dimensional reduction 
projects this equation onto another PDE, which governs 
the density p(x, t) in the reduced space of the coordinate 
x. The coordinate x is one of the coordinates of the full 
space, r = (x, y), and the mapping accomplishes integra- 
tion over the transverse coordinates y. In the case of the 
FP equation (jl.2[) . phase space (x,v) represents the full 
space and the dimensional reduction integrates over the 
"transverse" coordinate, the velocity v. Hence we have 
the defining relation between the densities p and p: 



p(x,t) 



p(x, v, i)dv 



(2.1) 
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The phase-space density p(x,v,t) is expected to be near 
the thermodynamic equilibrium and so approximately 
proportional to exp(— /3mv 2 /2), which provides conver- 
gence of the integral in Eq. (|2.1j) . 

Then the first step of the mapping is also integration 
of Eq. (|1.2p over v. If completed, we get 



(d t + vd x ) p(x, v, t)dv = 0; 



(2.2) 



the other terms are zero due to p(x, v, t) — > in the limit 
v — > ±oo. This equation represents nothing but mass 
conservation in the reduced space, 



d t p(x,t) + d x j(x,t) = 0, 



(2.3) 



where the ID flux j is defined by the relation 

j(x,t) = / vp(x,v,t)dv . (2.4) 



In contrast to diffusion, where j = —Dod x p is fixed, 
here the flux j is a function formally independent of the 
density p, so we also need the evolution equation for this 
quantity. After integration of Eq. (|1 .2[) multiplied by v, 
we obtain 

f°° U'(x) 
d t j(x, t)+d x / v 2 p(x, v, t)dv + —^-p(x, t) = 



/3m 2 



- 0mv2/2 d v e 0mv2/2 p{x,v,t)dv (2.5) 



[after some algebra and applying the definitions (|2.1I) and 
(12.41). This step recalls the Grad's method of moments 
[18lll9|. For the ID Kramers equation (|1.2[) . taking only 
a couple of the zeroth (p) and the first (j) order moment 
of the phase space density p is satisfactory for generat- 
ing the selfconsistcnt system of the mapped (real space) 
equations jU]) and ([27g|l.] 

The next key point of the mapping algorithm is that 
of expressing the full-space density p(x, v, t) using the ID 
density p(x,t) and also the flux j(x,t) in this case. This 
relation enables us to complete integrations in Eq. (I2.5[) 
and get the evolution equation for j(x,t), together with 
the ID mass conservation, in closed form. The initial 
task is to find the zeroth order approximation, valid in 
the limit m — > 0. Our first proposal for such a relation 
between p{x, v, t) and the reduced space quantities p and 
j is based on the following physical construction: 

For an infinitely small mass of the particle, the stochas- 
tic force thermalizes its velocity v almost immediately 
after any move along the spatial coordinate x. Similar to 
the transverse equilibration of the 2D density of a parti- 
cle diffusing in a narrow channel with biasing transverse 
force 15, 20- 22|, one could try the formula with sepa- 
rated Boltzmann factor in the fast relaxing "direction" 
v, p(x,v,t) ~ y/(3m/2np(x, t) exp(— (3mv 2 /2). It is easy 
to check that it does not work here; the flux j becomes 



zero according to Eq. (|2.4p . To prevent this failure, let 
us suppose that the distribution in v is shifted by the 
local mean (macroscopic) velocity vq, depending on the 
local flux, j(x,t) = vo(x,t)p(x,t). Then we have 



I f3m 
2^T 



[1 + j3mvvQ 



-/3mv 2 /2 



p(x,t). 



Retaining only these two terms in the square brackets 
and replacing vqp by the flux j, one gets 

p (x, V, t) = ^ e -0mv-/2 t) + pmvj ^ f)] ; 

(2.6) 

which will be taken as the sought zeroth order relation 
between p and the reduced space quantities, p and j. 

This heuristic formula will be verified later by the ex- 
act mapping algorithm. Still, one can check immediately 
that the relation (|2.6p represents correctly a kind of back- 
ward mapping of the ID (spatial) functions p and j onto 
the phase space densities p; if substituted for p(x,v,t) in 
the defining relations (|2.1j) and (|2.4j) . we obtain identi- 
ties. Applying Eq. (|2.6|) to Eq. (|2.5|l . the integrals over 
v can be completed and the result, 

d t j(x,t) +d?^- + ^-p(x,t) = -lj(x,t), (2.7) 



together with the mass conservation, Eq. (|2.3|) . forms a 
closed couple of PDE, governing the mapped quantities 
p and j. 

In the limit m — > 0, the first term in Eq. (I2.7|) . dtj, 
becomes negligible and the equation expresses the zeroth 
order relation between the flux j and the density p, 



1 



(2.8) 



If combined with the mass conservation, Eq. ()2.3|) , we get 
the Smoluchowski equation (11.31) ; l/f3-f = Dq represents 
the diffusion constant. 

The calculation presented shows that the Smolu- 
chowski equation (jl.3|) is related to the FP equation (jl.2|) 
in the same way as the Fick- Jacobs equation ll| to the 
diffusion equation valid in a narrow 2D channel. Both 
mapped equations describe an asymptotic behavior of the 
full-space density infinitely rapidly equilibrating in the 
transverse direction; the velocity v plays the role of the 
transverse coordinate for the FP equation. Our consider- 
ations indicate that the mass of the particle, m, becomes 
the small parameter, controlling the series of corrections 
to the Smoluchowski equation in the case when the trans- 
verse equilibration is not infinitely fast. 



Following the scheme of the mapping procedure [13 . 
14] , the next point is that of searching for the true re- 
lation between the phase space density p(x,v,t) and the 
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ID quantities p(x,t) and j(x,t), replacing the heuristic 
formula (|2.6p . valid for nonzero m. Without losing gen- 
erality, it can be written in the form 



Cj{x, v)p(x, t) 
+/3mvfi(x,v)j(x,t) . (2.9) 



If the operators Cj and f] (with d x implicit) are expandable 
in m, one can substitute for p(x, v,t) in the FP equation 
(|1.2p and fix the coefficients of these operators to satisfy 
this equation in each order of m, similar to the map- 
ping of diffusion. Then, using the relation of backward 
mapping (|2 .9[) in Eq. (|2.5p gives the expansion of the 
evolution equation for j and finally, in combination with 
mass conservation (j2.3p . the sought series of corrections 
to the Smoluchowski equation (11.31) in terms of the finite 
mass to. 

To verify whether this scheme is viable, we analyze the 
exactly solvable case with U(x) = in the next Section. 



III. EXACTLY SOLVABLE MODEL 



The exact solution of the FP equation (jl.2l) with no 
potential, U(x) = 0, is presented in this Section. We 
demonstrate the mapping on the example of the phase 
space density p[x, v,t) evolving from the initial density 
p(x, v, 0) with thermalized velocity v. The mapped equa- 
tion, as well as the form of the operators Co and f] in Eq. 
(|2.9[) can be found explicitly in this case. This analysis 
will direct us in construction of the recurrence mapping 
scheme in Sect. IV. 

The case U(x) — is exactly solvable [23j], the Green's 
function G(x, v, t; x' , v', t') of the FP equation p.2[) , 

dt+vd x - -L-d v e-^ mv2 / 2 d v e^ mv2 / 2 ] G(x,v,t;x',v',t') 
= S(x - x')S(v - v')S(t - t') , (3.1) 

can be derived explicitly (see Appendix A), 

7/3 6(r - t') 



G = 



in 



x exp 



r' - tanh(T - t')^/1 - q 2 
[2(£ - £') - (u + «') tanh(r - r')f 



4[t 



tanh(r — r')] 



[q(u 2 + u' 2 ) - 2uu'] 



1-q 2 

if expressed in the scaled coordinates, 

r = jt/2m, 
t = (/?m/2) 3/2 ^ 



,(3.2) 



(5m 2 



and q = exp[— 2(r — r')]. If the thermalized particle is 
inserted at time t = with a spatial distribution po(x), 



p(x,v,0) = J^p (x)e-e™ 2 / 2 , (3.4) 

V ^7T 

evolution of the density p is given by the formula 



p(x,v,t) 



j3m 
~2T 



dv' / dx'G(x,v,t;x',v',Q) 



xp (x')e^ W2/2 

Po (x>)dx> ( k - ef 



47rD r 



exp 



fl - e~ 2T ) 



2Q 



(3.5) 



Do = 1/7/3, the abbreviations 
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Q = t - - (1 - e- 
Z = Q - i (1 - e 



(3.6) 



are used and the integration over x' runs over the whole 
(unspecified) ID spatial domain. 

Then the spatial (ID) density p and the flux j are 
integrated directly according to Eqs. (|2.1I) and (|2.4I) . 



2 f P (x')dx' /Q 



xg-^e-K-Wfl, (3.7) 



= v, 



(3.3) 



p(x,t) 



It is easy to check that the mass conservation (|2.3p . 9 r p+ 
a/ /3m / 29^ j = in the scaled coordinates, is satisfied. 
The quantity Q plays the role of a "stretched" time |24| . 
For short times, t -C 2m/j, Q ~ r 2 and Z ~ 4t 3 /3. 
The formulas p.5p and (|3 . T[) describe correctly behavior 
of the Newtonian particles in this limit. The mapping 
procedure, as outlined in the previous Section, requires 
us to study asymptotic behavior in the opposite limit, 
t > 2m/7. 

For large times, r — > 00, the stretched time Q becomes 
r and the formulas <|3 . Tj) represent the general solution of 
the diffusion equation, as expected according to the pre- 
vious Section. Now it is necessary to verify that the mass 
m can serve as a small parameter controlling the series 
of corrections to the diffusion equation and its solution. 

It may look problematic at first glance, because Q 
contains cxp(— 2t) = exp(— "ft/m), representing essen- 
tial singularity of the variable m — > 0. Then the formu- 
las (|3.7p [and similarly Eq. (|3.5[) ] are not expandable in 
to. Nevertheless, this property is still consistent with the 
general scheme of the mapping, as analyzed in Ref. [3] ■ 

The dimensional reduction, as demonstrated on 
anisotropic diffusion in a narrow channel [3], also re- 
duces the Hilbert space of the full-space problem. Let us 
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denote M(e) the spatial operator of the evolution equa- 
tion, i.e. M(e) = 9^ + (l/e)<9^ for anisotropic 2D diffu- 
sion; the eigenvalues A.; and the eigenfunctions ipi(x,y) 
are given by the equation 



M(e)ipi(e;x,y) = \i(e)ipi(e;x,y), 



supplemented by proper boundary conditions at the walls 
of the channel. The parameter of anisotropy e < 1 splits 
the spectrum into two parts, the low-lying states, whose 
eigenvalues A;(e) remain finite for e — > and the tran- 
sients with A r (e) diverging ~ 1/e. Then the exact 2D 
density p evolves as 



P(e; x, y,t)=^r Cii>i(e; x, y)t 



-\i(e)t 



(3.9) 



the constants cj are given by the initial condition and 
the summation runs over the whole spectrum. The tran- 
sients contribute to the sum by the terms proportional 
to exp(— X r t/e), where A r = e\ r (e) are finite in the limit 
e — » 0. The result is a formula containing the essential 
singularity in the parameter e near zero, similar to Eq. 
p.5[) with singular exp(— jt/m) for m — > 0. 

The mapping procedure reduces the full Hilbert space 
of all ipi onto the space defined only by the low-lying 
states ipi. If the ID density p(x, t) is integrated from Eq. 
(I3.9[) and mapped backward onto the full Hilbert space 
(by some operator u>), the transients will be canceled; the 
sum (|3.9p after the mapping there and back runs only 
over the low- lying states ipi. The terms retained involve 
no essential singularity in e; the formula for p considered 
in the mapping procedure represents the regular part of 
the exact 2D density p with respect to the parameter e 
near zero. This reduction is natural for the zero-th order 
(Fick- Jacobs) approximation, as the transients decay in- 
finitely fast due to their infinite eigenvalues A r (e — > 0). 
Nevertheless, the mapping based on fixing the series of 
corrections expanded in e can work only with the regular 
part of the 2D density. 

Correspondingly, the formulas (|3.5p and (I3.7[) are ex- 
act, including the contributions of the transients, which 
are represented by the singular terms ~ exp(— yt/m). 
The mapping requires us to analyze only the regular 
parts, 



Preg t) 
jreg ^) 

and 



p(x')dx' 
2y/7rD (t - D Pm) 
x exp [-{x - x') 2 )/4L> (i - D j3m)] 
(x — x')p(x')dx' 



4^/irD (t-D l3m)3 
x exp [-(a; - x') 2 /AD (t - D f3mj] (3.10) 



p reg (x,v,t) = \ — 



(3m 



p(x')dx' 



x exp 



(x - x' - D^mvf 1 2 

— ■ — D 771 17 

AD (t -3D Pm/2) 2 H 



(3.11) 



written in the unsealed coordinates, obtained after taking 
only the regular parts of Q and Z (|3.6[) , Q reg = t — 



(3.8) 1/2 and Z reg = r - 3/4, in Eqs. (J3J]) and Of 
course, the formulas are applied for t 3> Doftm, when the 
transients vanish. We omit writing the subscript "reg" 
in the following text. 

In comparison with the overdamped limit m — > 0, evo- 
lution of the spatial (ID) density p and the flux j is only 
corrected by a time shift, t — > t — D f3m in the formulas 
(|3.10j) . The Gaussian distribution is retarded by the time 
to = D$fim = m/7, corresponding to the mean time nec- 
essary for losing information about the original velocity. 
The value of the shift to is constant in x and t, so the 
density p (|3.10|) still obeys the diffusion equation, 



d t p(x,t) = D d%p(x,t), 



(3.12) 



and j(x, t) = —Dod x p{x, t). There is no correction of the 
zero-th order mapped equation due to the finite mass of 
the particle. Let us stress that the case of U(x) = is 
extremely simple, similar to the mapping of the 2D dif- 
fusion in a flat narrow channel, also giving no corrections 
to the Fick-Jacobs approximation. 

Nevertheless, the relation of the backward mapping, 
generating the phase-space density p from the mapped 
quantities p and j, is not quite trivial. In the formula 
for p, Eq. (|3.11[) , the time is shifted by 3£o/2 and the 
displacement x—x' by vto with respect to the distribution 
of a massless particle. The shortest way to construct the 
relation of backward mapping is by applying the shift 
operators in t and x onp(x, t), compensating the different 
shifts of time and displacement in Eqs. (|3.10[) and (|3.1ip . 

pf x , v ,t) = x [^e- pmv2 ' 2 e~ {t °/ 2 ^~ vt ° d *p{x,t). 
V 2n 

(3.13) 

Due to Eq. (|3.12p . the operator dt in the exponent can 
be replaced by Dq8 2 . After expanding the shift operator 
exp(— vtod x ) and using j = —Dod x p, we arrive at the 
relation of the form (12.91). 



p(x,v,t) = \ —e 
Zir 



P m „-f3mv 2 /2-(D t a /2)dl 



x(D t ) k dl k 



k=0 

1 (3mv 



(2k) 



(2k + 1)1 



(3.14) 



the explicit formulas for Cj and r) to be substituted in Eq. 
(|2.9|) in the case U(x) = are 



2tt J 2y/nD (t-3D Pm/2) 



oj(x, v) 



fj(x,v) 



'1^ ( 2 fe)! {Doto) dx ' 



00 / q 2 \ k 

- (DoW2)a2 E|^T(^o) fe 9f.(3.15) 
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Both operators are expandable in m, (to = TO/7); their 
zero-th order coefficients equal unity, consistent with the 
heuristic formula (|2 .6[) . Also applying the relation (13. 14|) 
in the definitions (|2.1j) and (|2.4[) gives identity. 

The final step is that of verifying the evolution equa- 
tion for j, Eq. (|2.5p . Two integrals are to be completed 
with p expressed by the backward mapping, Eq. (|3.14p . 



v 2 p(x, y, t)dv = (l + D t Q dl) p(x, t)//3m, 

Xj 

^ rnv2 /2 d v e f}rnv2 12 p{x,v,t)dv = (3mj(x,t), (3.16) 



which are now valid exactly. If substituted in Eq. (|2.5I) . 
and the equation (|3. 12|) is applied, we get 



(l + t dt)\j{x,t) + D d x p(x,t)] =0 



(3.17) 



after simple algebra. This equation validates the relation 
j = —DodxP for nonzero m as well and thus, if combined 
with the mass conservation (|2.3|) , also the diffusion equa- 
tion (|3.12|) without any corrections. If compared with 
the calculation of Eq. (|2.8[) . the 1-st order term ~ d t j, 
neglected in the previous Section, is compensated here 
by other 1-st order term coming from the exact relation 
(|3.14j) . appearing in the integrals (|3.16[) . 

The FP equation (|1.2[) with zero potential is too sim- 
ple to give nonzero corrections to the diffusion equa- 
tion (|3.12[) . Nevertheless, it helped us to understand the 
structure of the mapping. It shows that the scheme sug- 
gested in the previous section is viable. The relation of 
the backward mapping has the form of Eq. (|2.9[) . at least 
for U (x) — 0, and the operators Cj and fj can be expanded 
in m, or to — TO/7, 



00 



lu = y^jfa k (x,u) 

fc=0 



fj = 



00 

= V" 



k—0 



(3.18) 



Integration in Eqs. (|3.16|) indicates that the coefficients 
and f/k should be sought dependent up on the scaled 
velocity u rather than v (compare to Ref. [IBj]); other- 
wise each term in Eq. (|3.18l) would contribute in several 
succeeding orders in the integrals (|3.16|) . The mixing 
of orders would hinder us in constructing the recurrence 
scheme generating corrections to Eq. (|2.5p . In the nota- 
tion of Eq. (|3.18[) , Cj — fjo — 1 and 



u)t(x,u) 
f]i(x,u) 

Cj 2 (x,u) 



(u 2 - l/2)D d 2 x , 
(u 2 /3-l/2)D d 2 x , 
(u*/6-u 2 /2 + l/8)D 2 di, 



(3.19) 



valid for U(x) — according to Eqs. (|3.15p . will be used 
for testing the results of the recurrence procedure in the 
next Section. 



IV. MAPPING PROCEDURE 

We now finish the construction of the mapping proce- 
dure, outlined in the Section II, for an arbitrary analytic 
potential U(x). Supposing the backward mapping of the 
form (|2.9p with the operators Cj and fj expanded in to 
(to) according to Eqs. (|3.18p . we find recurrence rela- 
tions fixing the coefficients Cj k and fjk- Completing the 
integrals in Eq. (|2.5p . we obtain a series of corrections to 
the zero-th order relation j = —Dod x p. Combined with 
mass conservation, Eq. (|2.3p . it gives the Smoluchowski 
equation corrected due to nonzero mass of the particle. 

The essential relation determining the operators Cj and 
fj is the FP equation (|1.2p . which has to be satisfied for 
any solution of the reduced problem, the density p(x, t) 
and the flux j(x,t), after their backward mapping (|2.9p 
onto the full-dimensional Hilbert space. If the expansion 
in to — to/7 of both operators, Eq. (|3.18p . is supposed, 
we have 



di 



2 Wg) 
■ — uo x a u 

/3m y/2PrH 



2tn 



-d u e~ u d u e u 



E*o 

fe=0 

xe"" 2 &k{x, u)p(x, t) + y / 2f3muijk(x, u)j(x, t)J =0(4.1) 

after introducing the scaled velocity u = y (3m /2v in Eq. 
(|1.2p . The factor /3m is replaced by to /Do in the following 
calculations. Thus half-integer powers of to appear in Eq. 
(|4.ip 25j ■ As this equation has to be satisfied for any to , 
we can split it for clarity into two relations: the first one, 
including only the integer powers of to, 



fc=0 



d t i2>kP + 2u 2 d x fjkj - f3U'(x)e u d u ue " fj k j 



zto 



-0, (4.2) 

and the second one, collecting the half-integer powers, 



2t ud t fikj + D (2ud x - (3U'(x)e u d u e 



xw fe p-e u d u e u d u ufj k j 



0. (4.3) 



Notice that Eqs. (14. 2p and (|4.3[) do not violate par- 
ity of cDfe and fjk in u. If used for construction of the 
recurrence relations between the coefficients, all they 
have to have the same parity as Cjo = fjo = 1; hence 
ojk(x,u) — (2>k(x,—u) and fjk(x,u) — fjk(x,—u). This 
symmetry enables us to find the normalization (or iden- 
tity) conditions for and fjk ■ The backward mapped p, 
Eq. (|2.9p . with the operators Cj, fj expanded in to, Eq. 
(|3.18jl . substituted in the definitions (flOj) and ([241) nas 
to give identities for any to, p(x,t) and j(x,t). Thus we 
obtain 



1 



due u Cjk(x, u) = 5o,k, 



(4.4) 



7 



2u 2 due u fj k (x,u) = S ,k- 



(4.5) 



The operators U)k and fjk are supposed not to depend 
on time, so the time derivative commutes with them and 
acts directly on p(x,t) or j(x,t) in Eqs. (|4.2j) . (|4.3I) . To 
derive the operators unambiguously, using only spatial 
derivatives, we express d t p = —d x j from the mass con- 
servation, Eq. (|2.3I) . However, the time derivative of j 
cannot be expressed in a similar way from Eq. (|2.5I) . 
because dtj is not the leading term there. If the back- 
ward mapping, Eqs. (|2.9|) and (|3.18j) . is applied to the 
integrals of Eq. (|2.5j) . we get 



v 2 p(x, v 



t)dv = ^~ f 1+ y>si fe (x) ]p(x,t), 



where the operators Ik are given by 
the right-hand side integral of Eq. (|2.5I) results in 



u 2 due u LOk(x, u); 



(4.6) 



7 



13m 2 
7 



■£*8 



due u d u ufik(x,u)j(x,t) 



j{x,t) 



to 



after integrating by parts and using the normalization 
relation (14.51) . Then, instead of the evolution equation 
for j, Eq. (|2.5[) has to be understood as an expression 
relating j and p, 



(l+t d t )j{x,t) 



-d x J2toh(x) 



fc=i 



p(a:,t), (4.7) 



and the flux j, as well as its time derivative, is expressed 
using p according to this equation. The term todt acts 
now on p as the 1-st order correction in (l+todt)~ 1 p after 
completing inversion and commutation with the spatial 
operators in the square brackets of Eq. (|4.7I) . Then dtp 
is systematically replaced by —d x j, contributing to the 
next corrections in the higher orders of to . The result is a 
formula for j expressed by some purely spatial operator 
acting on p, 



l + Y,toZk(x) 

fc=i 



d x eP u Mp(x,t), 



(4.8) 

where the operators Zfc are related to /fc, Eq. (|4.6I) . Using 
this form, we are able to write the final mapped equation 
explicitly after combination with Eq. 



d t p = D d x e 



-0U(x) 



1 



fc=i 



d x e^% 



^2^Q k p = Qp, 



(4.9) 



fe=0 



which is the Smoluchowski equation (jl.3l) in the zero-th 
order, extended by a series of mass dependent correc- 
tions in to. The operators Qk are introduced to simplify 
notation in the following calculations. 

Now the operators Zk can be expressed explicitly using 
Ik. Expanding (l + todt)~ 1 p in to and applying Eq. (|4.9|) . 
we have 



(l + t d t )- 1 p=(l-t Q + ...)p 



1 



k=0 



which is to be used in the operator equation, 



tnZk 



fc=l 



OC 



fc=l 



E*o +1 

fc=0 



QkY\ (4.10) 



obtained after comparison of Eqs. (|4.7p and (|4.8I) . Ex- 
panding in to an d comparing the coefficients of the same 
powers of to, we get a sequence of relations fixing Zk, 



e-^ZAeW 



dji - e^ u d x e' 3U Qo, 
d x l2 — d x I\Qo 
+ e-^d x e^{Ql-Q 1 ) 



(4.11) 



Finally, we derive the recurrence scheme, generating 
the operators uik, determining Ik, Eq. (14. 6[) . and thus 
also ifc or Qk, Eqs. (|4~9|) . (|4TTj) . Notice that the fj k do 
not directly enter the mapped equation (|4.9p . but they 
are necessary in the recurrence formulas for uik- 

The recurrence scheme for Cjk is defined by Eq. (|4.2p . 
This equation has to be satisfied for any p and j, solv- 
ing the mapped problem, but these quantities, although 
considered before as formally independent, are related by 
Eq. ()4.8|) . To obtain an equation for operators, j has to 
be expressed by p. Applying the relation (|4.8|) and Eq. 
(|4.9|) for d t p, we get 

OC 

J2 to D (u> n d x - 2u 2 d x f, n + pu'e^duue-^^y-^ 



xJ^toZAe^-^due- 

1=0 



1 

2to 



d u 0Jr, 



p = 0, (4.12) 



valid for any function p(x,t); we have Zq = 1. To low- 
est order, t^ , only the term e 11 9 tl e~ u "9„t2;o = 0. It 
is satisfied by wo = 1, the only solution nondiverging at 
u — > ±oo and also satisfying the normalization, Eq. ()4.4[) . 
In the higher orders, we derive the recurrence relation, 

n 

e u2 d u e' u2 d u Cj n+1 = 2L> £ ~~ 2y2d ^k + PU'(x) 

k=Q 
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xe u2 d u ue- u2 fj k )e-^Z n . k d x e^ x \ (4.13) 

Calculation of the uj n +\ requires us to know the fj k up 
to k = n. They are generated from Eq. f|4.3j) . Again, 
j, as well as dtj, have to be expressed by p using the 
relations P~g|) and (|43|) . Then Eq. flU} becomes 



n=0 L fc=0 

oo oo 

-2D o t o ufj n e- 0U 4Z k d x e 0u d x e- 0u Y] t Q Z x d x e &v 

k=0 



1=0 



+ (2ud x ~pU'e u2 d u e- u2 )to n 



p = 



(4.14) 



valid for any function p(x,t). In lowest order, t 



1/2 



e u2 d u e~ u2 duU^e-W d x e? u = -2ud x + (3U'e u2 d u e~ u2 
= -2u(d x + 0U') = -2ue- pu d x e^ u , (4.15) 

where we have used loq — 1 and Zq = 1. After the first 
integration, we have 



d u ufio 



2udue- u = 1 + Cie^ 



(4.16) 



the integration constant C\ — provides convergence as 
w — > ±oo. The next integration gives f/o = 1 + (1/u)Cq', 
Co = 0. This calculation validates our heuristic formula 
(|2.6p in the zero-th order approximation. 
In the higher orders, Eq. (|4.14[) generates the relations 

e u2 d u e- u2 d u ufj n e-P u d x eP u = ((3U'e u2 d u e~ u2 - 2ud a )u n 

k+Kn 

-2D Y, ufjn-k-t-te-^Zkd^d^-^Zid^ 



k,l=0 



n-1 

£ 

k=0 



e u ' 2 'd u e~ u2 'd u ufj k e- pv ' Z n - k d x e^ , 



(4.17) 



forming the recurrence scheme for r)„. Completing the 
operations in Eqs. (|4.17p and (|4.13p one has to keep 
in mind that the equation acts on an arbitrary function 
p(x,t), not depending on u. On the other hand, the 
operators Lu k and f] k for k > depend on u. 

The recurrence procedure starts from Cjq = 1 and Zq = 
1. Calculation of the next order correction requires first 
expressing the fj n according to Eq. (|4. 17[) . or (|4.15l) for 
ri = 0, as shown above. Then Cj n +i is derived from Eq. 
(|4.13[) . I n +i integrated according to Eq. (|4.6I) and finally 
Z n+ i expressed from Eq. (|4.11l) . To demonstrate the 
procedure, we derive the first order correction, Z\. 

We use already calculated r?o = 1. For n — 0, Eq. 
(|4.13[) becomes 

e u2 d u e~ u2 d u Qj x = 2£> (1 - 2u 2 )e~ Siu ^ ^^l. (4.18) 
After the first integration, 

dud = 2D e u2 (ue- u2 + d) e^^e^, (4.19) 



the integration constant G\ = 0, to provide convergence 
for u — > ±oo. The integration constant Co after the 
next integration is fixed to satisfy the normalization, Eq. 



due- u2 D Q (u 2 + CoK^^e^^ = 0, (4.20) 



hence Co = —1/2 and 

£>! = D a (u 2 - l/2)e-^dy u ^. 



(4.21) 



For U(x) — 0, we recover the corresponding formula in 
Eq. (|3.19p . Integration over u in Eq. (|4.6p results in 
h = Dq exp[— fiU(x)]d x exp[f3U(x)], giving finally 

Z x = D {e^d x e~^d x - d x e^ u d x e^ u ) = D (3U"(x) 

(4.22) 

fromEq. (|4~TT|) . 

In the higher orders, fj n , tl) n +\ arc calculated according 
to Eqs. (I4.17P and (|4.13l) . The integration constants 
after double integration have to provide convergence for 
u — > ±oo, requiring the operators to be even in u, and 
the normalization, Eq. (|4.4[) . The condition (|4. 5[) for 
fj n is satisfied automatically; it serves as a check on the 
computation. The derivation is tedious, and we present 
only the results in second order, 



= D (u 2 /3^1/2)e-^dy u ^\ 



6j 2 



Zo = =^ 



fU 



iH 3 



+ {u 2 - -) (4PU"(x)8 x + 3(3U^(x))\d x e^\ 

°° r 4(/3C/"(x)) 2 - p 2 U'(x)U^(x) + pu^ix) 
+3/3U {3 \x)d x ]. (4.23) 



Again, the formulas (|3.19l) for U (x) = are recovered. 
There are no contributions to Z n in this case, too, as ex- 
pected according to the analysis in the previous Section. 
Also, linear potentials, U(x) — —Fx, have no effect on 
validity of the uncorrected Smoluchowski equation (|1.3p . 
The particle driven by a constant force F move asymp- 
totically with constant mean velocity vq — F/j and the 
distribution p(y,t) in coordinate y, shifted by the drift, 
y = x — Vot, is again Gaussian as in the case of no poten- 
tial. 

The situation becomes different if the driving force 
F(x) is not constant. If the mass m or the time of the 
thermalization to — m/j is small, but nonzero, the par- 
ticle appearing at a new position x has to accommodate 
to the new local mean velocity. On the other hand, it 
carries some mean momentum from its previous position 
and needs some time to change it. Meanwhile it slips to 
some other position than predicted by purely stochastic 
dynamics due to its inertia, or non-zero mass. The ef- 
fects of such slipping are indicated by the corrections Z n 
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of the Smoluchowski equation and they are nonzero for 
potentials with nonzero U"(x), or higher derivatives. 

As seen from Eq. (|4.23[) , the Z n are not only functions, 
but operators, containing d x in the higher orders. So the 
mapped equation (|4.9I) has exactly the same structure 
as the mapped equations for diffusion [l6|, or biased dif- 



xiapp ec 

fusion 15J, |22|, |26|. Being inspired by these works, Eq. 



(I4.9[) can be simplified by replacing the correction oper- 
ators A)[l + J2kLi txjZk] by a function D(x), a spatially 
dependent effective diffusion coefficient, 

d tP (x, t) = d x e-^D(x)d x e^p(x, t), (4.24) 

which becomes valid in the limit of stationary flow, i.e. 
when the spatial density p and the flux j change very 
slowly, p(x,t) — > p(x). Due to mass conservation, Eq. 
(j2.3[) , the flux j(x,t) — j is constant (but nonzero) in x 
as well. If expressed from Eq. (|4.24[) . 



_ _ p -PU{x) 



D(x)d x 



*p(x), 



(4.25) 



the function d x (exp[(3U(x)]p(x)) — —jexp[f3U(x)}/D(x) 
is dependent only on the system; m, 7 and the poten- 
tial U(x), for any stationary solution p{x). So we can 
substitute for it in Eq. (|4.8|) . 



fc=i 



,f)u{x)_3_ 
D(xY 



(4.26) 



valid for stationary flow, and calculate D(x) unambigu- 
ously from the expansion of the corrections 



D(x) 
as a series in to 
D(x)/D 



fc=i 



(4.27) 



Dot Q pU" + (D Q t Q ) 2 [2((3U") 2 
+P 2 U'U {3} + l3U (4) /2\ + ... .(4.28) 



The next simplification is that of neglecting all the 
derivatives but U"(x), i.e. approximating the real U{x) 
locally by a quadratic potential. In this case, the expan- 
sion (|4.28l) can be summed up to infinity, 



D(x)/D 



(2n)! 



n\(n 

n=0 y 



1) 



[D o to0U"{x)] n 



ADoto(3U"(x) 
2Dot /3U"(x) 



(4.29) 



the proof is given in the Appendix B. The formula works 
for D t /3U"{x) = mU"{x)/ry 2 < 1/4, which is the con- 
dition for non-oscillatory movement of a particle in a 
quadratic well with friction, the damped harmonic os- 
cillator. If U"(x) — k is constant, the trajectory of a 
particle averaged over the stochastic force is governed by 



m(x) + j(x) + k(x) = 



(4.30) 



DM 

Do 



1.5 



-0.5 



-0.25 



0.25 



mU"(x)// 



FIG. 1: Plot of the effective diffusion coefficient D(x) depen- 
dent on the mass of the particle m and friction 7 according 
to Eq. (14.291) . valid if the derivatives of the potential higher 
than U"(x) are neglected. 



from Eq. p. II) . The particular solutions are (x(t)) = 
exp(ai) with a = —(jzt \/"f 2 — Amu) /2m. Requiring a 
to be a real number gives the same condition. 

This simple example demonstrates restriction of the 
theory presented to nonoscillatory movement of the par- 
ticle in potential wells on its way along a ID channel. A 
small mass m is expected, to enable the friction quickly 
to damp the momentum of a particle; i.e. to have the 
relaxation in the velocity faster than in the real space 
coordinate x. 

A more detailed insight to the restrictions of the di- 
mensional reduction of the phase space controlled by the 
mass m can be obtained by comparison of the Green's 
function (GF) of the mapped equation (|4.24j) with D(x) 
given by Eq. (|4.29p and GF of the Kramers equa- 
tion (jl.2p for the damped harmonic oscillator, U(x) — 
kx 2 /2 = mLu 2 x 2 /2, which is exactly solvable. The solu- 
tion G = G{x, v, t; x' , v' , t') of the equation 



(d t + vd x - uj 2 xd v 



7 



/3m" 



■d v 



~pmv 2 /2Q e l3mv 2 /2^jQ 



S(x - x')S(v - v')S(t - f) (4.31) 



reads [l[ 



G 



■ exp 



(si - s 2 )e^/ m 
2nVab - h 2 

+2h(£ ~ £o)(v - Vo))/2(ab - h 2 )] , (4.32) 



where 



Sl,2 



2m 



± 



Am 2 



n 

a - 



(s\x — v)e S2t 
(s2X — v)e~ Slt 

— ~- 2s " 



fim 2 s\ 



£0 
Vo 
b = 



h 



2 / 
/3m V 



= (six - 

~- (S2X 1 - 
1 

fim 2 S2 

„7*A 



(4.33) 



-2s 2 t 



(4.34) 
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Similar to the case of no potential in the Section III, 
let us suppose that a thermalized particle (equilibrated 
in velocity) was inserted at a position Xq at time t' = 0, 



p (x\v') = ^6(x> -x )e-P™' 2 / 2 . (4 



35) 



After integration over v' and v, we get the corresponding 
spatial density 



p(x, t) = I G(x,v,t;x'v' ,0)po(x' ,v')dx'dv'dv 



(3m 
2^~Z 



wo(si - s 2 )exp 



2Z 



Z = 



:( Sl (x - x e S2t ) - s 2 (x - x e stt fj ; (4.36) 



. Sl (l + e S2t )- S2 (l + e Slt ) Sl (l- e S2t )- S2 (l-e Slt ) 



On the other hand, the coefficient D(x), Eq. (|4.29l) . 
becomes constant for the quadratic potential, 



D(x) 



"f/2m — \J "f 2 /Am 2 — coq 



si 



(4.37) 



/3mu>Q fimujQ ' 

and GF of the corresponding mapped equation (|4.24p . 

(d t - d x D{x)e- f)m ^ x2 l 2 d x e fSm < x2 l 2 )g{x,t ] x M) 

= d(x-x )S(t-t ), (4.38) 

can be easily found by a calculation similar to the deriva- 
tion of Eq. (|3.1[) . Appendix A. The result, 



g(x,t;x ,0) 



/3m 



2tt(1 



3 2 Slt ) w oexp 



1 



-(3muj 



x(a;-a; e Sl *)7(l-e 2sit ) , (4.39) 



describes evolution of the real space density of a ther- 
malized particle inserted at xq, too, and can be directly 
compared with the formula (|4.36l) . 

First, let us notice that in comparison with Eq. (|4.36|) . 
the exponential e S2t disappeared from the formula (|4.39|) . 
The root s 2 — — 7/m for m — > makes e S2t ~ e -7 */" 1 the 
term essentially singular in m and so invisible for the 
recurrence procedure, which works with the operators Cj 
and f] expanded in m. Using our argumentation from 
the Section III, e S2 * represents the "transients" neglected 
by the mapping. On the other hand, e Slt ~ e -™"o*/7 [ s 
regular in m small, representing the contribution of the 
low-lying states, retained by the method. 

Next, let us stress that the equation (|4.24[) with D(x) 
expressed by the expansion (|4.28|) was derived in the limit 
of the stationary flow; for the net flux almost constant, 
which is not the case of the process described by the 



Eqs. (|4.36[) and (|4.39|) . Nevertheless, using the approxi- 
mations 



si 



s 2 (l± 



( Sl -s 2 )(l±e s 



S\X - 



s 2 (x-x e Slt ) ~ (si-s 2 )(a 



x e 



applicable for e Sl * -C 1, the regularized formula (I4.36P 
(with e S2t neglected) becomes finally Eq. (|4.39|) : i.e. it 
represents correctly the asymptotic behavior of the spa- 
tial density p(x, t) for large time t. 

If m approaches 7 /2uio, the transients contributing by 
e S2 *, neglected by the mapping, become important. In 
the oscillating regime, si j2 are complex numbers and both 
are necessary for expressing the real density p(x, t) in 
Eq. (|4.36|) . The mapping which splits the Hilbert space 
to the retained low-lying states and the neglected tran- 
sients, controlled by m small, loses its justification and 
the method stops working. Mapping in this region re- 
quires a different method to be applied. It will be an 
object of our study in the future. 



V. CONCLUSION 

Although modeling of transport in confined systems is 
often based on study of the Langevin equation, Eq. (|1.1[) 
in the simplest ID case, the solutions necessary in prac- 
tical applications are often accessible only in two limits: 
either the friction 7—^0, when the particles obey New- 
tonian dynamics, or the mass m — > 0, which corresponds 
to stochastic dynamics. Any solution in the region be- 
tween these limits requires working in phase space, which 
makes the problem much more complicated. 

The present paper shows how to describe the region 
of finite 771/7 while still working in real space, as in the 
case of stochastic dynamics. The equation governing evo- 
lution of the spatial density p(x, t) is the Smoluchowski 
equation, corresponding to the limit of a massless par- 
ticle, extended by a series of corrections in powers of 
to = m/7; to can be interpreted as the typical time 
of thermalizing of the particle's initial velocity by the 
stochastic force. 

In general, the extended Smoluchowski equation has 
the form of Eq. (|4.9|) . Do = l/jf3 denotes the diffusion 
constant and the operators Zk are systematically derived 
within the recurrence procedure presented in the Sect. 
IV. In the limit of stationary flow, i.e. when the flux 
is almost constant but nonzero in time and space, this 
equation can be simplified to Eq. (|4.24|) . where the ef- 
fective diffusion coefficient D(x) is calculated unambigu- 
ously from the operators Zk, Eq. ()4.28|) . In the simplest 
approximation, when all the derivatives of the potential 
higher than U"(x) are neglected, the series of corrections 
can be summed up to infinity, giving the formula for D(x) 
in a closed form, Eq. (|4.29|) . described in Fig. 1. Then 
the equation describes stationary flow in a quadratic po- 
tential. The theory works while AmU"{x) < 7 2 , until the 
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averaged trajectory of a single particle is not oscillatory 
in the potential wells along the ID channel. The map- 
ping in the oscillatory regime requires the next analysis, 
which will be done in the future. 

Technically, the paper demonstrates that the projec- 
tion technique developed for mapping of diffusion in 2D 



can 



(3D) channels with varying cross section [13. Il4 
be adapted for the dimensional reduction of a process de- 
scribed by an evolution equation of a different type than 
the diffusion or Smoluchowski equation. The method has 
been modified significantly; m/7 had to be confirmed as 
the small parameter controlling the expansion of the cor- 
rections Z kl as well as the operators of the backward 
mapping, Cj and r). In contrast to diffusion, the flux 
j{x,t) is handled here as a quantity independent of the 
density p(x,t). Thus the recurrence procedure, calculat- 
ing expansions of the correction operators Z k and the 
operators a), r), is in principle the result of combination 
of three equations, Eqs. (jl.2[> . (|2.3I) and also (|2.5|) . with 
the relation of the backward mapping, Eq. (|2.9[) . 

It is worthwhile to notice that including the mass de- 
pendent corrections to the Smoluchowski equation re- 
sults in the equations (|4.9[) or (|4.24l) . which are of the 
same form as the comparable equations obtained from 
the mapping of diffusion in channels with varying cross 
section. On the other hand, the effective coefficient D(x) 
(I4.29[) has a different symmetry than the similar formu- 
las extending the Fick- Jacobs equation 1(| 17 1 for con- 
fined diffusion. Also D(x) can be greater than 1 here 
(see Fig. 1); i.e. the quasi stationary flux is accelerated 
when passing through a shallow potential well, depend- 
ing on the nonzero mass of the particles. These interest- 
ing properties could be observed in simulations similar to 
that verifying D{x) in the extended Fick- Jacobs equation 



27 



28j . The effects of slipping of the particles diffusing 
under a nonconstant force F{x) due to their inertia, as 
described in Sect. IV, might also influence the interesting 
phenomena in the micro and nano world, such as Brow- 
"29I 30j . rectification of the flux in quasi ID 
stochastic resonance or the 



man pumps 
structures 



31] 



negative mobility [35|. Study of such applications of the 
theory presented is expected in the future. 
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APPENDIX A: EXACT SOLUTION 



The Green's function (|3.2j) solving the FP equation 
with zero potential, Eq. (|3.1[) . is calculated here. First 
we introduce the scaled coordinates £, u, r according to 



Eqs. p. 31) and define the function r(£, u, r; u',t'), 

G(x, v, t; x', v', t') = e-" 2 / 2 r(£, u, r; u', r')e u ' 2/2 , 

(Al) 

satisfying the transformed equation (|3. 1[) . 
[d T + udt - d 2 u + u 2 - 1] r(£, «, r; r') 

= ^e^-^'HiZ - £')<*(« - u')6(t - r'),(A2) 

the exponential factor becomes unity due to 6(u — vf). 
After the Fourier transform in £ and r, 

(A3) 

and shifting the velocities by ik/2, w = u + ik/2 and 
w' = u' + ik/2, the equation 

[-iv -d 2 w +w 2 -l + k 2 /4] T ktU (w; w') = ^-6(u - u') 

(A4) 

becomes solvable if Tk^(w;w') is expressed in the basis 
set of the linear harmonic oscillator ip n (w), 

00 

) = ^r n (M)V>nMCK)- (A5) 

n=0 

The eigenfunctions ip n (w) satisfy 

(-dl + w 2 ) 4> n (w) = \ n ip n {w) = (2n + l)i(t») (A6) 

and we use the integral representation of the Hermite 
polynomials H n {w) [36| 



^/7rV2" 



\=H n (w)e' w2 ' 2 



2" 



Itt 3 / 2 



tiItt' 



e- w2 / 2 / (w + ir) n e- r2 dr (A7) 



in the next calculation. 

Using the transformations above, we find 



T n (k,v) = 



7/3/4 



< + A„ - l + fc 2 /4' 



(A8) 



Applying it in the formulas (|A5|) and (|A3|) . we integrate 
the last one over v in the complex plane, 



SttDo 

OO 

x^e-^-^MwWniw'), (A9) 

n=0 

S(x) denotes the Heaviside unit step function and Dq = 
1/7/3 is the diffusion constant. Now the integral relation 
(|A7|) is used for i/) n (w) and ipn( w ') ancl the summation 
over n can be readily completed. Finally, the straightfor- 
ward triple integration over k, r, r' is performed and us- 
ing the transformation (|A1|) results in the formula 
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APPENDIX B: QUADRATIC APPROXIMATION 



Derivation of the formula (|4.29[) for the effective dif- 
fusion coefficient D(x) with all the derivatives higher 
than U" (x) neglected is presented here. This approxi- 
mation corresponds to local replacing of the potential by 
a parabola, U(x) ~ k(x — xq) 2 /2 + Uq, where re, xq and 
Uq are fitting parameters. 

First we simplify Eq. (|4.27[) . For quadratic potential, 
the right hand side can be rewritten as 

e-^W (l + toZ) ~* e' 3U ^ = (l + t e-P u Z( X y u ) , 

t()Z = Y^n=i toZ n ', the difference contains only the higher 
derivatives of U(x), which are zero. Hence 



D(x)/D = l + e-( }U WJ2 t oZne 



PU(x) 



(Bl) 



n=l 



The formulas for D{x) have been derived considering 
stationary flow; j(x,t) — j is constant. It simplifies the 
relation (|4.7|) ; dtj = 0. Thus the right hand side rep- 
resents stationary flux, which can be directly compared 
with Eq. (|4.8p . giving a much simpler relation between 
Z n and I n than Eq. (|jTTD|) . 

e-P u ^Z n (x)d x eP u ^p(x) = dj n (x)p(x) (B2) 

for any stationary solution p(x). Calculation of the coef- 
ficients of D(x) according to Eq. (IB1[) requires us to take 
d x exp[/3U(x)]p(x) = exp[/3U(x)], hence finally 



corresponding to the normalization of f) n , Eq. (|4.5[) , and 
the relations p,16[) . p. 171) . proving no correction to the 
Smoluchowski equation in the case U (x) = 0. 

The operators ui n and fj n for the quadratic potential 
have the form 



k=l 
n 



Vn 



^ U{X) E Cn,kQ k (u) {pU'T k dfe^\ (B6) 
fc=i 

with the coefficients 

2fc (2ti-1)! 



c n,fe — f o 



(n - k)l(n + k)\ ' 



(B7) 



Due to the integrals, Eq. (|B5j) . only the first terms with 
Pi (it) in Eq. (IB6|) contribute to the expansion of D(x), 
Eq. (|B3I) . Then the functions become 



e -pu{x)z e pu(x) 



2 (2n- 1)! 



(n - l)!(n+l)! 



(DqPU'T, (B8) 



taking J7^ 3 ^(x) = into account. Applied in Eq. (|B1[) it 
results in the expansion of D(x), Eq. (I4.29[) . 

Finally, one has to verify that the formulas (|B6[) sat- 
isfy the recurrence relations (14. 13)) and (|4.17l) . acting on 
the function p(x) — exp[— /3U(x)] J dx exp[(3U (x)} . Al- 
though the equations simplify notably due to neglecting 
the derivatives higher than U"(x), we omit the details of 
this tedious but straightforward calculation. 



V 77 J-oo 

XLu n (x,u)e- 0u( - x) f dxe pu{x) (B3) 



after application of Eq. (|4.6[) . 

Before writing the explicit formulas for Qj n for the 
quadratic potential, we define the polynomials 

n I i \n- fco2fc — n 



w«) = E 



fc=0 



(2fc)!(n- fc)! 

(-1) 
(2fc + l)!(n- fc)!' 



n — k<r\2k — 7 



2k 



(B4) 



n = 1, 2, coming from the expansions of w„ and r) n for 
zero potential in to, Eqs. p. 151) . The first few polynomi- 
als are visible in the round brackets of Eq. p.!9[) . One 
can check by direct integration that 



Q n (u)u 2 e u du 



i+T l^i 



(-If 



k—r 



2 n+i ^ kUn-k)\ 

k=0 y ' 



= 0, 



P n {u)u e u du = —S n ,i, 



(B5) 
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